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ABSTRACT 


Microwave  interactions  with  inhomogeneous  plasmas 

are  often  studied  "by  employing  a  simplified  electromagnetic 

approach,,  i.e.,  by  representing  the  effects  of  the  plasma 

by  an  effective  dielectric  coefficient.  The  problems  and 

approximations  associated  with  this  procedxjxe  will  be 

discussed.  The  equation  describing  the  microwave  field 

in  an  inhomogeneous  partially  ionized  plasma  will  be 

derived,  and  the  method  that  has  been  applied  to  obtain 

the  reflected,  transmitted,  and  absorbed  intensities  in 

inhomogeneous  plasmas  will  be  presented.  The  interactions 

of  microwaves  with  plasmas  which  have  Gaussian  electron 
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density  profiles  have  been  considered  by  Klein. 

In  his  work,  as  well  as  in  other  papers  treating  micro - 
wave  interactions  with  inhomogeneous  plasmas,  the  variation 
of  collision  frequency  with  position  has  been  neglected. 

In  general,  the  assumption  of  constant  collision  frequency 
is  not  justified;  e.g.,  for  a  highly  ionized  plasma, 
the  electron  density  profile  determines,  in  part,  the 
profile  of  the  electron- ion  collision  frequency.  The 
effect  of  the  variation  of  the  collision  frequency 
profile  on  the  interaction  of  microwaves  with  inhomo¬ 
geneous  plasmas  has  been  studied  in  order  to  obtain  an 
estimate  of  the  degree  of  error  which  may  resuJLt  when 
constant  collision  frequency  is  assumed  instead  of  a 
more  realistic  collision  frequency  profile.  It  will  be 
shown  that  the  degree  of  error  is  of  particular  importance 
when  microwave  analysis  is  used  as  a  plasma  diagnostic. 
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MICROWAVE  IMERACTIONS  WITH  INHOMOGENEOUS  PARTIALLY  IONIZED  PLASMAS 
A.  H.  Kritz,  Space  Science  Laboratory,  General  Dynamic s/Astronautics 

I.  INTRODUCTION 

Theoretical  considerations  indicate  that  if  the  nature  of  a  plasma 
is  sufficiently  simple^  microwave  interactions  with  the  plasma  may  yield 
diagnostic  information.  The  interactions  of  microwaves  with  the  charged 
plasma  particles  are  normally  treated  hy  describing  the  plasma  in  terms 
of  an  effective  dielectric  coefficient.  The  assimiptions  made  in  formula¬ 
ting  the  dielectric  treatment  will  he  outlined  below  in  Section  II.  In 
Section  III^  the  analyses  and  approximations  employed  in  interpreting 
the  microwave  data  will  be  reviewed.  Also,  the  results  of  some  of  the 
previous  theoretical  studies  of  inhomogeneous  plasmas  will  be  summarized. 

The  effect  of  the  collision  frequency  profile  on  microwave  plasma 
diagnostics  will  be  discussed  in  Section  IV.  In  the  final  section,  the 
problems  associated  with  utilizing  the  diagnostic  theory  for  microwave 
interactions  with  inhomogeneous  partially  ionized  plasmas  will  be 
summarized.  Also,  the  possibility  of  using  electromagnetic  scattering 
from  electron  density  fluctuations  as  a  possible  diagnostic  tool  will  be 
discussed  briefly. 

II.  THE  EFFECTIVE  DIELECTRIC  COEFFICIENT 

In  the  derivation  of  the  effective  dielectric  coefficient  for  a  partially 
ionized  plasma,  the  following  assumptions  are  required:  The  electrons, 
ions,  and  neutrals  are  taken  as  three  independent  gases  with  interactions 
smoothed  into  continuous  forces.  Only  friction-like  action  and  reaction 
(collisional)  forces  and  macroscopic  electromagnetic  forces  are  considered. 
Mean  velocities  are  employed  in  writing  the  equations  of  motion  for  the 
three  gases,  and  non-linearities  are  omitted  by  a  perturbation  treatment. 

Since  the  gas  must  be,  to  a  close  approximation,  electrically  neutral, 
the  average  ion  density  is  taken  equal  to  the  average  electron  density 
(n^  =  ^  simple  dielectric  treatment  would  not  be  possible  without 

the  additional  assumptions  that  the  velocity  of  the  plasma  as  a  whole 
is  zero  and  that  the  gradients  of  the  electron  and  ion  pressures  are  zero. 
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In  order  to  "be  able  to  solve  the  equations  of  motion  for  the  three 
gases^  it  is  necessary  to  obtain  the  expressions  for  the  friction-like 
collision  forces.  The  collisional  force  exerted  on  the  j *th  particle 
gas  by  the  k’th  particle  gas  is  given  by  the  average  momentum  a  typical 
j  P'crticle  loses  per  collision  times  the  number  of  collisions  per  second, 
V.,  .  The  j-k  particle  collision  frequency  is  given  by 


y  =  Q  I V .  -  V  '  n,  =  (n, /n.)'^^.  • 

j  k  ^ j  k  1  ^k  k  ^  k^  j  ^  k j 


(1) 


where  Q..  is  the  cross  section  for  the  j-k  particle  collision:  v.,  n., 

and  m.  are  the  velocity,  nmber  density,  and  mass,  respectively,  of 
J 

■type  j  gas.  The  subscripts  e,  i,  and  a  will  be  used  to  denote  the  electron, 
ion,  and  neutral  gases,  respectively.  Since  m^  «  m^,  the  electrons  lose, 
on  the  average,  a  quantity  of  momentum  equal  to  their  mean  momentm 
Therefore,  the  rate  of  loss  of  momentum  per  unit  volume  by  the  electron  gas 
due  to  collisions  with  the  ions  is 


F  . 

'--ei 


-m  n  V  . ( 
e  e  ei 


(2) 


The  rate  at  which  the  electron  gas  loses  momentum  is  equal  to  the  rate 
at  which  the  ion  gas  gains  momentum  from  the  electron-ion  collisions, 
i.e.,^^^  =  Similarly  for  electron-neutral  collisions. 
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v)==nv  m(v  -  v)=-F 


(3) 


The  positive  ions,  in  collision  with  the  neutral  atoms  of  the  same  mass 

(it  is  assumed  M  =  m.  =  m  ),  on  the  average,  lose  one-half  their  momentum 

X  a 

relative  to  the  neutral  gas.  Therefore, 
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On  the  basis  of  elementary  theory,  the  following  equations  of  motion 
for  the  electron,  ion,  and  neutral  gases  can  now  be  written  employing 
the  expressions  for  the  frictional  forces  developed  above: 
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n  e 
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n  m  V  =  n  eE  +  — ^  (v  xH)+F.+F 
e  e-^e  e  ^  c  we  ^ei  -^ea 

#  ^  ^ 

n  Mv.  =  -n  eE  -  —  (v.  x  H  )  +  F,  +  F, 
e  e  ^  c  Vi  ^le  -^la 


n  Mv  =  F  +  F  . 
a  -'a  ^ae  '^ai 


(5) 

(6) 

(7) 


where  E  is  the  applied  microwave  fields  and  is  a  constant 
field  Imposed  on  the  plasma.  The  microwave  fields  are  taken 
harmonic  time  dependence  e*"^^  so  that  =-icj^j. 


magnetic 
to  have  a 


Equations  (5)  to  (7)  can  now  be  solved  for  the  electron  and  ion 
velocities,  and  v^,  in  terms  of  the  applied  microwave  field, 
consequently,  the  complex  conductivity  tensor,  ^p,  relating  the  current 
density,  j,  and  the  electric  field  can  be  obtained: 


j  =  en  (v  “  V. )  =  a  •  E. 
e^^e  1 


(8) 


Consider  now  Maxwell's  equations; 

=  ikH 


VxE  = 


1 

c 


,  „  *  1  -  *  “t* . 


-ikf I  +  a) 

O)  ' 


(9) 


=  -ike^  •  ^ 


(10) 


where  k  =  cd/c. 
e,  is  given  by 


From  Equation  (lO)  it  is  seen  that  the  dielectric  tensor. 


■t-ni 

O) 


s) 


(11) 


where  o  is  obtained  from  Equation  (8).  (The  magnetic  permeability  has 
been  taken  as  unity  in  Maxwell's  equations  since  the  permeability  has 
been  shown^  to  differ  only  infinitesimally  from  unity  for  most  plasmas 
of  physical  interest).  If  the  external  magnetic  field,  H^,  is  assumed 
to  vanish,  the  dielectric  tensor  can  be  treated  as  a  scalar,  and  Equation 
(lO)  can  be  written 


V  X  H  =  -ikeE 


(12) 


If  the  motion  of  the  ions  is  also  neglected,  then  the  dielectric 
coefficient  has  a  fairly  simple  representation, 

(O)  /od)^  (v/cd)(cd /cd)^ 

e  =  1 - P - ^  +  1  - E - ^  (13) 

1  +  (v/o))  1  +  (v/a>) 

where  v  =  v  .  +  v  and  o)  is  the  plasma 
ei  ea  p 
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It  is  seen  here  that  under  the  approximations  outlined  above  the 
difference  in  the  treatment  of  a  partially  ionized  plasma  as  to  a  fully 
ionized  plasma  is  reflected  in  the  value  of  the  collision  frequency; 
i.e.,  in  a  fully  ionized  plasma  the  electron-neutral  collision  frequency 
vanishes,  =0,  In  the  inhomogeneous  plasma  both  o)^  and  v  are 
functions  of  position.  Equations  (9),  (12),  and  (13 )  provide  the  basis 
for  the  approximate  electromagnetic  treatment  of  a  partially  ionized 
inhomogeneous  plasma, 

II.  INTERACTIONS  WITH  INHOMOGENEOUS  PLASMA5 

The  following  equations  for  the  microwave  electric  and  magnetic 
fields  are  obtained  from  Equations  (9)  and  (12): 


frequency  given  by 
1/2 


(iJ+) 


+  k'^eE  =  - 

(15) 

2  2 
.  H  +  k  eH  =  - 

(  2^  )  X  (  V  X  H). 

(16) 

If  the  plasma  is  uniform  (  V  e  =0),  Equations  (I5)  and  (I6)  reduce  to 


+  k^eE  =  0 

(17) 

+  k^€H  =  0 

(18) 

where  e  is  given  by  Equation  (13)  and  is  a  function  of  position  in  an 
inhomogeneous  plasma.  Also,  if  the  variation  in  electron  density  and 
collision  frequency  is  normal  to  the  electric  field  (E  •  V  e  =  O),  the 
electric  field  satisfies  Equation  (l?),  or  if  the  variation  is  in  the 
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direction  of  E  (E  x  V  €  =0),  the  magnetic  field  satisfies  Eq^uation  (18). 

Equations  (17)  and  (I8)  describe  the  microwave  fields  within  the 
plasma  when  the  geometric  optics  approximation  is  made^  i.e., 


“  I  «  k.  (19) 

Under  this  approximation,  the  right  hand  sides  of  Equations  (15)  and  (16) 
are  of  small  order  of  magnitude  so  that  they  can  be  neglected.  For  most 
applications  of  geometric  optics  €  >  1;  consequently,  the  condition  expressed 
in  Equation  (I9)  is  satisfied  when  V  e  is  comparatively  small.  However, 
when  the  plasma  is  described  in  terms  of  an  effective  dielectric  coeffi¬ 
cient,  €  may  be  less  than  unity,  and  in  the  region  e  ^  0^  the  condition 
expressed  in  Equation  (19)  can  not  be  satisfied  even  when  V  €  is  small. 

An  alternate  procedure  for  the  solution  of  the  microwave -plasma  interaction 
problem,  based  on  the  integral  representation  of  the  field  equations  and 
the  application  of  the  Born  approximation,  will  be  presented  below.  First, 
however,  some  theoretical  studies,  based  on  the  differential  field  equations, 
will  be  discussed.  In  these  studies,  semi-infinite  plasmas  with  electron 
density  variations  along  the  rectangular  coordinate  normal  to  the  plane 
plasma  boundary  are  considered;  thus,  for  normal  incidence,  ^  *  7  e  =  0, 

2 

V.  A.  Bailey  has  given  a  summary  of  the  approximate  analytic  tech¬ 
niques  that  have  been  applied  to  solve  an  equation  having  the  form  of 
Equation  (I7)  (with  e  a  function  of  position).  In  addition,  Bailey 
presented  a  new  approximate  solution  which,  although  somewhat  cumbersome, 
avoids  some  of  the  difficulties  to  which  the  more  customary  approximate 
solutions  are  subject.  Also,  the  solutions  obtained  by  Bailey's  technique 
appear  to  have  a  greater  range  of  validity  than  previous  approximate  solu- 
t ions . 

There  have  been  several  studies  of  inhomogeneous  plasmas  in  which  the 
electron  density  is  chosen  so  that  the  solution  to  the  field  eqxiatlons 
could  be  expressed  in  tenns  of  common  special  functions.  Pappert  and  Plato 
considered  a  linear  electron  profile  (the  solution  in  terms  of  Hankel 
functions  of  one-third  order),  a  parabolic  profile  (the  solution  in  terms 
of  Hermite  functions),  and  a  cosine  profile,  i.e.,  the  argument  varying 
from  zero  to  n/2  (Mathieu  functions).  They  discussed  the  results  of  the 
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linear  profile  in  some  detail  especially  with  application  to  the  reentry 
telemetry  problem.  The  reflection  and  transmission  of  electromagnetic 
waves  at  linear  electron  density  gradients  has  been  studied  also  by 
Albini  and  Jahn.  Taylor  shoved  that  the  solution  to  the  field  equation 
could  be  expressed  in  terms  of  cylindrical  Bessel  functions  for  both  an 
exponential  rise  in  electron  density  and  an  electron  density  varying  as 
X  .  Although  certain  electron  density  profiles  can  be  treated  analytically, 
the  results  of  studies  of  the  interactions  of  microwaves  with  inhomogeneous 
plasmas  having  special  electron  density  profiles,  such  as  those  listed 
above,  have  limited  utility  since  (l)  only  very  few  profiles  can  be  treated 
analytically,  (2)  the  reflection  and  transmission  coefficients  contain 
rather  complex  functions  which,  in  general,  are  not  tabulated,  and  (3) 
the  variation  of  collision  frequency  with  position  is  neglected  (either 
constant  collision  frequency  is  assumed  or  the  collision  frequency  is 
omitted  entirely). 

In  a  paper  by  Klein,  Greyber,  King,  and  Brueckner^  and  in  several 
7-10 

reports  by  Klein  a  series  of  problems  involving  the  interactions  of 
microwaves  with  inhomogeneous  plasmas  are  treated  by  employing  the  WKB 
approximation  and/or  numerical  methods;  reflection,  transmission,  and 
absorption  amplitudes  are  presented  for  many  cases.  In  the  p ape r^  and  in 
one  of  the  reports,  an  exponential  electron  density  profile  is  considered; 

n(x)  =  n  exp  ^-(x/d)  x  >  0, 

where  d  is  measured  in  units  of  free  space  wavelength  and  the  microwaves 
are  Incident  from  the  left.  In  an  attempt  to  imderstand  the  interaction 
of  microwave  radiation  from  a  re-entry  vehicle  with  the  stirrounding  plasma,  ^ 
non-normal  incidence  is  treated  by  employing  geometrical  ray  tracing  proce¬ 
dures  which  restricted  the  calculations  to  short  wavelengths.  In  another 
report, microwaves  are  assumed  incident  from  the  right  on  an  exponential 
profile  (valid  for  all  x).  In  order  to  investigate  the  interaction  of 
microwave  radiation  with  a  non-uniform  plasma  sheath  adjacent  to  a 

9 

conducting  wall,  Klein  treated  an  exponential  profile  and  a  parabolic 
profile  adjacent  to  a'  conductor.  The  Gaussian  electron  density  profile, 

r  J 

n(x)  =  exp  -  (x/d)  j  ,  (20) 
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is  considered  in  detail  for  the  purpose  of  understanding  microwave  inter- 
actions  with  a  non-uniform  ionized  wake.  In  this  study,  as  in  the  studies 
described  above,  constant  collision  frequency  is  assumed.  The  Gaussian 
density  profile  was  employed  in  our  investigation  of  the  effect  of  the  vari¬ 
ation  of  the  collision  frequency  on  the  interactions  of  microwaves  with 
inhomogeneous  plasmas  which  is  presented  in  Section  IV.  Klein  investigated 
the  possibility  of  replacing  the  Gaussian  electron  variation  by  a  homogeneous 

slab  of  thickness  2d  and  electron  density  equal  to  n  .  He  found,  however, 

max 

that  the  rectangular  electron  barrier  was  a  good  approximation  only  if  the 
slab  is  very  thin  compared  to  the  free  space  wavelength,  i.e,,  d  <  .05. 

When  the  slab  thickness  is  of  the  order,  or  greater  than,  the  free  space 
wavelength,  order  of  magnitude  errors  in  the  reflection,  transmission,  and 
absorption  amplitudes  may  result. 

Maxwell's  equations  in  differential  form  (Eqs.  I5  and  I6)  have  also 

been  applied  in  the  study  of  the  reflection  of  radio  waves  from  the 

ionosphere.  In  these  studies  it  is  necessary  to  include  the  effect  of 

the  earth's  magnetic  field.  Comprehensive  reviews  of  this  problem  have 
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been  given  by  Budden  and  Ginzburg. 

Because  of  the  complexity  of  the  problem  of  microwave  interactions 
with  inhomogeneous  plasmas,  when  plasmas  involve  coordinates  other  than 
rectangular,  it  has  been  generally  necessary  to  resort  to  expressing 
Maxwell's  equations  in  integral  fomi  and  to  applying  the  Born  approxima¬ 
tion.  The  interaction  problem  for  a  few  electron  profiles  (such  as 
radial  variation  proportional  to  r^  for  cylindrical  and  spherical 
coordinates)  can  be  solved  analytically  using  the  differential  form 
of  Maxwell's  equations.  However,  more  often  than  not,  the  solutions 
involve  special  functions  which  are  not  tabulated.  Below,  we  shall 
discuss  briefly  the  integral  representation  of  the  field  equations  and 
some  studies  based  thereon. 

The  differential  field  equations,  Eqs.  (I5)  and  (I6),  can  be  written 
in  the  following  form: 

V  X  V  X  E  -  k^E  =  k^(e-l)  E  (2l) 

V  X  V  X  H  -  k^H  =  -Ik  V  X  (€-1)^  ,  (22) 
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By  considering  the  right  side  of  Eq.uations  (21)  and  (22)  as  source  terms, 
the  standard  Green’s  function  technique  can  be  applied.  Using  the  free 


space  Green's  function, 


exp  [^lfc(r,-r ' )][ 


(23) 


where  ^  and  r’  are  the  coordinates  of  the  source  and  observation  points, 
respectively,  we  can  express  Eqs .  (2l)  and  (22)  in  integral  form: 

2^,  ,  .n 


E(r)  =  E^(r)  +  jjj  (r  ’  )-l]  f  ^  E(r  ’ )  •  V’]  V'G(r,  r  '  )+k  E(rJ  )G(r,r  ’  )j  dv  ’ 


H(r)  =  H^(r)-iic  V  X  JJJ  ' 


(24) 

(25) 


where  E  and  H  are  total  fields  and  E  and  H  are  the  incident  fields. 
If  the  Herztian  potential  due  to  induced  polarization  [€(rj-lj  E(r) 
is  expressed. 


[€(r')-3  E(^)G(r,r')dv', 


(26) 


then  Maxwell’s  equations  can  be  written  as  follows; 


E  (r)  =  E-E  =  V  V  •  n  +  k St 


(27) 


H  (r)  =  H-H  =  -IkV  x  n  .  (28) 

—  s  -^~o  - 

When  the  condition  ^€(r)-l^  «  1  is  valid  for  the  plasma  under  considera¬ 
tion,  the  Born  approximation  can  be  applied  to  Equations  (24)  and  (25), 
i.e.,  ^(ij)  in  the  integrand  can  be  replaced  by  E^(^').  This  implies 
that  the  field  to  which  the  scattering  electrons  are  subject  is  only 

the  incident  field.  In  general,  for  this  condition  to  be  valid,  the  electron 

2  2 

density  has  to  be  very  low  (w  «  u)"")  or  the  collision  frequency  very  high 
2  ^ 

(vo)  »  0)  ),  i.e.,  absorption  is  neglected;  consequently,  the  range  of  validity 
of  this  approximation  is  quite  limited. 
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The  integral  formulation  was  used  by  Bart  he  1  to  calculate 

scattering  from  a  narrow  cylindrical  column,  a  homogeneous  cylinder. 
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and  a  cylinder  with  a  Gaussian  radial  distribution.  Barthel  has  also 
considered  the  inverse  problem^  i.e.,  the  determination  of  the  electron 
density  distribution  when  the  scattered  field  is  known  as  a  function  of 
k  =  Chu  and  Politis^^  in  calculating  the  radar  cross  section  of 

the  ionized  wake  of  a  hypersonic  body  used  the  field  equations  in  integral 
fonii.  They  considered  a  cylindrical  plasma  having  an  electron  density 
distribution  of  the  form 

n.  =  n  exp  ]  -  — rr  -  ~  .  ( 29 ) 

max  ^  b  -  ^ 

With  this  density  distribution  and  the  application  of  the  Born  approxima- 
tion_,  the  integrations  in  Eq.  (2^)  can  be  carried  out  in  closed  form. 
Results  are  presented  for  various  experimental  conditions. 

Another  procedure  for  obtaining  diagnostic  information  from  the 
interaction  of  microwaves  with  inhomogeneous  plasmas  will  be  discussed 
in  Section  V.  First,  however,  we  shall  present  some  results  from  a  study 
of  the  effect  of  the  variation  of  collision  frequency  on  microwave  inter¬ 
actions  with  inhomogeneous  plasmas. 

IV.  EFFECT  OF  THE  COLLISION  FREQUENCY  PROFILE 

The  studies  of  inhomogeneous  plasmas  have  indicated  that  if  micro - 
wave  analysis  is  to  be  used  as  a  plasma  diagnostic,  the  correct  electron 
density  profile  must  be  used  in  the  calculations.  However,  in  the 
previous  studies,  constant  collision  frequency,  though  not  justified, 
is  assumed.  In  some  plasmas  (e.g.,  highly  ionized  plasmas  where  the 
electron  density  determines,  in  part,  the  collision  frequency  profile) 
the  collision  frequency  profile  may  have  an  important  effect  on  the  micro- 
wave  interactions.  Because  the  variation  of  collision  frequency  with 
position  has  been  neglected  and  because  only  a  limited  number  of  density 
profiles  have  been  treated,  a  study  was  undertaken  to  develop  an  efficient 
computer  program  for  determining  the  microwave  interactions  with  semi¬ 
infinite  plasma  slabs  which  have  arbitrary  electron  density  and  collision 
frequency  profiles.  The  program  has  been  applied  to  plasma  slabs  which 
have  Gaussian  electron  density  distributions  in  order  to  investigate  the 
effect  of  varying  the  sharpness  of  Gaussian  shape  collision  frequency 
profiles . 
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The  width  of  the  semi-infinite  inhomogeneous  plasma  slab  is  determined 
by  the  criteria  that  outside  the  slab  the  magnitude  of  the  plasma  dielectric 
coefficient  differs  from  the  free  space  coefficient  (i.e.^  unity  for  the 
Gaussian  units  used)  by  less  than  one -one  thousandth.  The  object  of  this 
study  has  been  to  obtain  an  estimate  of  the  degree  of  error  which  may 
result  when  constant  collision  frequency  is  assumed  instead  of  a  more 
realistic  collision  frequency  profile.  Knowledge  of  the  degree  of  error 
is  of  particular  importance  when  microwave  analysis  is  to  be  used  as  a 
plasma  diagnostic,  for  it  will  be  shown  that  the  collision  frequency  pro¬ 
file  may  have  a  rather  dramatic  effect  on  the  interactions  of  microwaves 
with  inhomogeneous  plasmas.  A  preliminary  report  of  these  results  has 
been  presented  elsewhere. 

The  plasma  electron  density  and  collision  frequency  profiles  which 
we  have  considered  are 


03  /go)  2 

=  r-(x/d)2] 

(30) 

v/o.) 

"  exp[-'Nx/d)^j  . 

(31) 

The  width  of  the  Gaussian  density  distribution  is  determined  by  d 
(measured  in  units  of  the  free  space  wavelength,  the  velocity  of  light 
divided  by  the  applied  frequency);  the  parameter  N,  which  we  vary 
between  zero  and  two,  determines  the  sharpness  of  the  collision  frequency 
profile.  The  range  of  the  parameters  which  we  considered  are  as  follows: 

0.l6  <  (OD  /cu)^  <  100 

—  p'  max  — 

.001  <  (v/o))  . 

.01  <  d  <  3 

The  reflection,  transmission,  and  absorption  of  the  microwaves  by  the 
plasma  slab  were  obtained  by  integrating  the  differential  field  equation 
(Eq.  17)  numerically.  The  technique  used  in  the  integration  is  a  backward 

17 

difference  method  suggested  by  Cowell.  The  difference  equation  is  based 
on  Stirling's  interpolation  formula  and  does  not  contain  the  first  or  any 
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odd  differences.  Consequently,  the  integration  procedure  is  rather 
eTficient,  A  brief  outline  of  the  integration  procedure  is  presented 
in  Appendix  I. 

An  amplitude  and  phase  for  the  wave  transmitted  by  the  slab  are 
assumed,  and  the  usual  boundary  conditions  of  continuous  tangential 
fields  are  employed.  The  computer  program  yields  the  incident  and 
reflected  waves  (and  associated  phase  angles)  corresponding  to  the  assumed 
transmitted  wave.  The  reflection  and  transmission  amjilitudes  are  then 
readil;y  obtained.  When  this  study  was  started,  it  was  hoped  that  rules 
of  thumb  could  be  developed  for  taking  into  account  the  collision 
frequency  profile  so  that  embers ome  nmerical  caJ.culations  could  be 
avoided.  This  objective  has  been  accomplished  only  in  part,  for,  within 
the  range  of  variables  considered,  the  collision  frequency  profile  often 
has  a  rather  dramatic  effect. 

The  microwave  interactions  with  plasma  slabs  having  four  different 
collision  frequency  profiles  will  be  presented  in  order  to  illustrate 
the  effect  of  the  variation  of  collision  frequency  on  the  interaction  of 
microwaves  with  the  inhomogeneous  plasmas.  The  four  profiles  are  shown 
in  Figure  1.  The  profiles  considered  are: 

(v/oi)  =  exp  j^-(Nx/d)^^ 

N  =  0  Constant  collision  frequency 

N  =  0,5  Broad  Gaussian  profile 

N  =  1.0  Gaussian  profile  equivalent  to  that  of  the 

electron  c^ensity 

N  =  2.0  Sharp  Gaussian  profile. 

In  Figure  2,  the  reflected  intensity  (normalized  to  \init  incident  in¬ 
tensity)  is  plotted  as  a  function  of  d  (one-half  the  distance,  in  vinits  of 
free  space  wavelength,  between  the  points  at  whicn  the  electron  density 
falls  off  to  one  e*th  its  peak  value).  In  the  plasma  referred  to  in 

Fig.  2^  ^  =1.0.  It  can  be  seen  that  significant 

p  max  max 

errors  in  plasma  diagnosis  would  result  if  variations  in  the  collision 
frequency  profiles  are  neglected.  In  Fig.  3^  "the  microwave  reflection 
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is  plotted  as  a  function  of  (od  /cjo)"^  for  a  plasma  characterized  by 

^  p'  ^max 

d  =  0.3  and  (v/cd)  =10,  Here  also,  it  is  apparent  that  diagnostic 
max 

calculations  of  electron  density  would  be  incorrect  if  constant  collision 
frequency  is  assumed  for  a  situation  when  the  collision  frequency  profile 
is,  in  fact,  Gaussian.  Note  that  the  N  =  0.5  profile  yields  a  rmaller 
reflected  intensity  than  when  the  collision  frequency  is  constant.  The 
reflected  intensity  first  decreased  as  N  varied  from  0  to  0.6  and,  then, 
the  reflected  Intensity  reversed  sharply,  and  increased,  as  N  continued 
to  become  larger.  The  cause  of  this  behavior  is  as  follows:  As  N  varied 
from  0  to  0,6,  the  Gaussian  profile  became  sharper,  lowering  the  average 
collision  frequency.  As  a  consequence,  the  reflection  from  the  first 
surface  increases,  but  the  rate  of  microwave  absorption  in  the  plasma  also 
increases  so  that  the  microwave  intensity  reflected  from  the  second  surface 
is  decreased  to  the  point  that  the  net  result  is  a  decrease  in  total  reflected 
intensity.  For  a  semi-infinite  plasma  (as  opposed  to  a  plasma  slab) 
decreasing  the  collision  frequency  produces  a  monotonic  increase  in 
reflected  intensity.  When  N  is  greater  than  0.6  practically  all  the 
reflection  occurs  at  the  first  sizrface  so  that  the  reflection  increases 
with  N. 

Attempts  were  made  to  obtain  criteria  for  deteznnining  effective 
constant  collision  frequencies  to  replace  the  Gaussian  profiles.  However, 
in  some  instances,  for  fixed  d  or  for  fixed  peak  plasma  frequency,  the 
reflected  intensity  is  increased  more  by  replacing  the  constant  collision 
frequency  by  the  sharp  Gaussian  (with  N  =  2)  than  by  replacing  it  with  a 
constant  profile  of  one-tenth  the  original  value;  whereas,  in  other  cases, 
the  constant  profile  of  one-tenth  the  original  collision  frequency  resulted 
in  the  larger  increase  of  the  reflected  intensity.  For  the  range  of  d 
considered,  it  was  also  noted  that  the  difference  in  the  reflected  intensity 
resulting  from  the  collision  frequency  profiles  corresponding  to  N  =  1 
and  N  =  2  decreases  and  becomes  insignificant  as  the  peak  plasma  frequency 
increases  and  the  peak  collision  frequency  decreases.  This  is  shown  in 
Table  1  below: 
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Insignificant  Difference  in  Reflected  Intensity 
from  N  =  1  and  N  =  2  Profiles  for 


(v/o)) 

"  max 

=  0.01 

and 

All  (CD  /o)' 
p  ■ 

2 

)  considered 

max 

( v/o)) 

^  '  "max 

=  0.1 

and 

(co  /o))^ 

^  p"  max 

>  10 

(v/cd) 

^  '  "max 

=  1.0 

and 

(cD  /cd)^ 

^  p"  max 

>  40 

(v/ce) 

^  "  "max 

=  10 

and 

(oD  /cu)^ 

^  p"  "max 

>  X  where  X  >  100 

Table  1 

Figures  k  and  5  demonstrate  the  effect  of  the  same  four  Gaussian 
collision  frequency  profiles  those  shown  in  Fig.  l)  on  the  trans¬ 

mission  of  microwaves  through  the  inhomogeneous  plasma.  In  Fig.  4  the 
transmitted  intensity  is  plotted  as  a  function  of  slab 

in  which  d  =  0.3.  In  general,  for  the  peak  value  of 

(cd^/od)  <1.0  and  d  <  1.0,  the  difference  in  the  transmitted  intensity 
resulting  from  a  constant  collision  frequency  profile  and  the  profile  for 
N  =  0.5  is  not  significant.  Figure  5  shows  the  transmission  as  a  function 

of  d  for  the  plasma  in  which  (o)  /od)  =  2.0  and  (v/o))  =  10.  From 

^  p'  ^max  max 

Figs.  4  and  5  it  is  seen  that  although  the  error  which  would  result  in 
the  transmitted  intensity,  if  an  incorrect  collision  frequency  profile 
were  assumed,  is  not  as  dramatic  as  the  resulting  errors  in  the  reflected 
intensity,  the  errors  in  the  calculated  transmitted  intensity  would  be 
significant.  The  following  two  conclusions  about  the  reflected  and  trans¬ 
mitted  intensities  can  be  made  from  this  study:  (l)  Within  the  range  of 
variables  considered,  when  (v/oj)  <  .01,  the  effect  of  the  Gaussian  collision 
frequency  profile  is  not  significant;  and  (2)  the  importance  of  using  the 
correct  collision  frequency  profile  increases  with  increasing  d. 

Finally,  we  shall  briefly  consider  the  effect  of  the  collision 

frequency  profile  on  microwave  absorption  in  an  inhomogeneous  plasma.  In 

Fig.  6,  the  absorbed  intensity  is  planned  as  a  function  of  d  for  a  plasm' 

2 

in  which  (od  /gd)  =  40  and  (v/oo)  -  1.0.  When  the  peak  value 
'  p'  "max  ^  "  "max 
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(cjd^/cjd)  >  1  the  effect  of  the  collision  frequency  profile  is  rather 
dramatic  as  shown;  whereas  for  ^  absorption  simply 

increases  with  d.  For  the  range  of  d  and  (cn  /<n)  considered,  when 
the  peak  value  of  (v/oa)  <  1  the  absorption  decreases  as  N  increases, 
i.e,,  as  the  collision  frequency  profile  becomes  sharper. 

In  investigating  the  importance  of  using  the  correct  collision  fre¬ 
quency,  a  computer  program  has  been  developed  for  determining  microwave 
interactions  with  plasma  slabs  which  have  arbitrary  electron  density  and 
collision  frequency  profiles  so  that  we  are  not  restricted  to  a  small 
nTJunber  of  profiles . 

V.  SUMMARY 

By  the  presentation  in  the  previous  section  of  the  sample  results 

indicating  the  effect  of  the  collision  frequency  profiles  on  plasma- 

microwave  interactions,  we  have  attempted  to  demonstrate  tliat  for  certain 

ranges  of  plasma  parameters,  the  collision  frequency  profile  may  have  a 

very  important  effect  on  the  interactions  of  microwaves  with  plasmas.  For 

example,  from  the  data  presented,  it  is  apparent  that  order  of  magnitude 

errors  in  the  electron  density  and/or  plasma  thickness  may  resiilt  if  the 

variation  of  collision  frequency  with  position  were  neglected.  From  Klein*s 

studies,  as  well  as  the  other  studies  of  inhomogeneous  plasmas,  it  can  be 

seen  that  except  when  either  the  electron  density  variations  occur  over 

distances  very  small  compared  to  the  wavelength  or  when  the  electron  density 

is  very  low,  incorrect  results  would  be  obtained  if  a  uniform  plasma  with 

some  average  electron  density  were  assmed  to  represent  the  inhomogeneous 

plasmas.  Thus,  since  a  variation  of  electron  density  or  collision  frequency 

with  position  may  yield  results  quite  different  than  those  obtained  when 

a  uniform  plasma  is  assumed,  extreme  care  must  be  exercised  when  micro- 

wave  analysis,  based  on  an  effective  dielectric  description,  is  applied 

for  plasma  diagnosis.  For  example,  the  anomalous  results  obtained  when  micro- 

l3 

waves  are  reflected  from  the  wake  of  a  re-entering  body  may  be  a  consequence 
of  neglecting  the  variation  of  electron  density  and  collision  frequency. 

V/e  have  not  yet  applied  the  resul.ts  of  our  study  to  practical  problems  for 
our  first  object  was  to  determine  the  importance  of  using  correct  profiles. 
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Before  concluding,  we  shall  briefly  mention  the  possible  utilization 

of  electromagnetic  scatter  from  density  fluctuations.  R.  Pappert  recently 

considered  the  application  of  this  technique  with  regard  to  analysis  of 

19 

the  wake  of  a  hypervelocity  projectile,  and  a  nimiber  of  groups  are 
currently  engaged  in  the  development  of  such  a  laboratory  tool.  The  the¬ 
ory  for  electromagnetic  scatter  has  been  considered  in  connection  with 

20-23 

ionospheric  studies.  When  the  effective  dielectric  theory,  described 

in  Sections  II  and  III,  is  employed,  it  is  possible  to  probe  the  ionosphere 

only  to  altitudes  corresponding  to  the  F  layer;  however,  by  observing  the 

backscatter  from  electron  density  fluctuations,  it  is  possible,  in  principle, 

24 

to  probe  beyond  the  peak  ionospheric  electron  density.  Recently,  Bowles 
reported  successful  probings  of  the  ionosphere  and  extensive  tests  are 
planned  for  the  future. 

The  technique  of  employing  electromagnetic  scattering  as  a  plasma 
diagnostic  tool  is  based  on  the  following  considerations:  Under  appropri¬ 
ate  conditions  (i.e.,  1,  the  scattering  volume  dimensions  large  compared 
with  the  incident  radiation  wavelength  but  small  compared  with  the  distance 
to  point  of  observation;  2,  the  Bom  approximation  valid;  3^  the  inter - 
electron  spacing  large  compared  with  scattered  wavelength;  and  4,  the  time 
of  observation  long  compared  with  ratio  of  scattered  wavelength  to  electron 
rms  thermal  speed),  it  can  be  shown  that  the  number  of  electrons  in  a  scatter¬ 
ing  volume  can  be  obtained  by  measuring  the  intensity  of  scattered  radiation, 
the  distance  from  the  scattering  volume  to  the  point  of  observation,  and 
the  angle  between  the  incident  and  scattered  radiation.  Also,  if  the 
primary  broadening  mechanism  is  Doppler  broadening,  the  scattered  Intensity 
per  unit  frequency  interval  may  yield  either  the  electron  or  the  ion 
temperatures.  If  the  inter-electron  spacing  is  small  compared  to  the 
scattered  wavelength  (the  opposite  of  the  limit  considered  above),  the 
differential  cross  section  per  unit  frequency  yields  the  Fourier  transform 
of  the  density  space  time  auto-correlation  function  from  which  it  may  be 
possible  to  infer  the  electron  density  distribution. 

I  wish  to  acknowledge  the  aid  given  by  Adolf  Hochstim  and  Arthur 
Anderson  in  the  investigation  of  the  effect  of  collision  frequency  profiles 
on  microwave  plasma  diagnostics.  Also,  I  wish  to  express  thanks  to  Richard 
Marriott  for  having  indicated  the  applicability  of  the  numerical  method 
used  in  the  calculations. 
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APPENDIX  I 


The  microwave  fields  propagating  in  an  inhomogeneous  plasma  in 
the  direction  of  the  gradient  of  the  dielectric  coefficient^  i.e.,  V  g(x)^ 
satisfies  the  equation 

^  e(x)E  =  0,  (l-l) 

dx 

Because  the  dielectric  coefficient  for  a  plasma  is  a  complex  function, 

e(x)  =  e^(x)  +  l€^(x)  (1-2) 

(see  Eq.  13 ihe  study  of  an  inhomogeneous  plasma  requires  the  simultane¬ 
ous  solution  of  the  following  two  differential  equations: 

p 

u"  =  (e  .u  -  G^v) 

1  i 

v"  =  -4jT^(e^u  +  (l-3) 

where  the  independent  variable  is  in  units  of  free  space  wavelength  and 
E  =  u  +  iv.  Equations  (1-3)  can  be  solved  by  a  straightforward  extension 
of  the  procedure  outlined  below  for  the  Integration  of  a  second  order 
differential  equation  of  the  form 

y"(T)  =  f(T)y(T)  +  g(T).  (1-4) 

The  solution  of  Eq.  (l-4)  is  based  on  the  following  difference  equation 
derived  from  Stirling's  interpolation  formula: 


^1 ' 


1 


+ 


(1-5) 


where  h  is  the  interval  by  which  t  is  increased  at  each  step  of  the  integra- 
tion  (taken  such  that  5  y"  and  higher  order  terms  are  negligible)  and 
y^  is  the  value  of  y  at  the  i*th  step.  The  difference  notation  for  a 
function  y  evaluated  at  constant  interval  can  be  understood  from  Table  I-'’ 
below: 
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.-1 


^1/2  ^\/2 


&  yi 


&  y. 


6“3y3/2  b-V3/2  &V3/2  &^y3/2 


0-2  .2 
B  72  &  72 


p,-l 

^  >^5/2  ^>^5/2 


Table  I-l 


where  the  terms  are  connected  by  the  relationships 


^n+1 

^  7. 


a/2 


"  ^  ^l/2(s+l)  "  ^  ^l/2(s-l) 


(1-6) 


for  any  integer  s.  By  substituting  y  from  Eq.  (1-4)  into  Eq.  (I-5)  and 
solving  for  y,  one  obtains 


=  h%  [6“^  y. 


i  12  ®il 


(1-7) 


a. 


=  1  -  f 

1  12  i  • 

Once  tvo  starting  values  for  y  are  known,  y^  and  y^>(ln  the  plasma  study 
where  only  one  value  of  y  is  known,  the  value  of  y,  a  distance  h  away,  is 
generated  by  a  Taylor’s  expansion)  the  nmerical  solution  to  Eq.  (I-3) 
can  be  tabulated  as  follows: 
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T 

a. 

1 

"1 

Si 

yV 

5.-2  tt 

&  y^ 

^■"■^1+1/ 2 

0 

a 

0 

f 

0 

So 

si  “'2  tt 

^  ^0 

(1) 

h 

“l 

Si 

^1 

yi' 

Si-2  II 

6  y^l^ 

(ii) 

2h 

“2 

"2 

S2 

(iv) 

(v) 

(Hi) 

(vi) 

3h 

“3 

■F 

3 

®3 

4h 

Slf 

Table  1-2 
^2 

The  known  functions  f,  g,  and  ^  tabxzlated  at  the  required 

interval.  The  values  for  y"  and  y"  can  be  obtained  by  substituting 

^  -P  ^  -P 

y^  and  y^  into  Eq.  1-4.  Then  6  y^’^  6  y”  are  calculated  from  Eq.  (l-?)* 

The  integration  routine  then  proceeds: 


(i) 

(ii) 

(iii) 

(iv) 


6'^  yi’/2  =  ^1'  ■ 


^3/2  =  yv2  ^ 


Vl  =  y’’  + 

Y?  r.-2  „  ^  1 

V  =  —  '  o  V  +  -  B  ' 

J'2  cCg  ^2  12 


(v)  yg  =  +  82 


(vi)  5'^  y'^I2  =  ^3/2  ^2 
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MICROWAVE  REFLECTION 


II 


Lr»  C3  o 
■  ■  ■ 

O  ^  CM 

II  II  II 
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2h 


MICROWAVE  TRANSMISSION 
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v/6)=(v/6>)„^  exp  [-(Nx/df 


II 
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MICROWAVE  ABSORPTION 
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(units  of  free  space  wavelength) 


